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Nomenclature 


A 

A* 

a"^ 

a“ 


A 


X 


a 


a(y) 

A 


B 

B, 


X 


c* 

V 


c* 

p 


D 

E 

F 


h(y) 


Helmholtz free energy. 

A/P V. 
c 

power law amplitude for specific heat above (see Table I). 

power law amplitude for specific heat below (see Table I) . 

parameter in eq. (3.13) for crossover function F. 

linear model equation of state parameter (see eq. (2.14c)). 

scaling function for A* (see Table II). 

coefficients in eq. (3.33) for ri of steam. 

o 

coefficients in eqs. (3.18) and (3.23) for of carbon 
dioxide and steam. 

power law amplitude for coexisting densities (see Table I) . 

parameter in eq. (3.13) for crossover function F. 

linear model equation of state constant (see eq. (2.14a)). 

coefficients in eq. (3.32) for r| of steam. 

coefficients in eq. (3.22) for X of steam. 

heat capacity at constant volume. 

C T /P V. 

VC c 

heat capacity at constant pressure. 

C T /P V. 
pc c 

specific heat at constant volume, 
specific heat at constant pressure. 

power law amplitude for chemical potential at T = T^ (see Table I). 

X/pc = thermal dlffusivlty. 

P 

constant in MLSG equation of state (see eq. (2.10)). 

crossover function in thermal conductivity equation 
defined by eq. (3.13). 

scaling function for y* (see Table II) . 


v 



= 

p ^(8p/9P)^ = compressibility. 


k 

= 

linear model equation of state parameter (see eq. (2. 

14b)). 


= 

Boltzmann's constant. 


n 

= 

exponent in eq. (3.13) for crossover function F. 


P 

= 

pressure. 


P 

c 

= 

critical pressure. 


P* 

= 

P/P . 
c 


P*(T*) 

0 

= 

analytic background in the equation for the pressure 
critical region (see Table II) . 

in the 

p*i 

= 

coefficients in eq. (2.17) for P*(T*). 


p(6) 

= 

auxiliary function in linear model equation of state 
(see Table III). 


Pi 

= 

coefficients in equation for p(9) (see Table III). 


q 

= 

parameter in eq. (3.27) for critical viscosity enhancement. 

R 

= 

universal constant in eq. (2.21) for 5 • 

o 


r 

= 

parametric variable in linear model equation of state 
(see eq. (2.14)) . 


S 

= 

entropy. 


S* 

= 

ST /P V. 
c c 


s(0) 

= 

auxiliary function in linear model equation of state 
(see Table III) . 


s'(e) 

= 

ds(0)/d0. 


S . 

1 

= 

coefficients in equation for s(0) (see Table III). 


T 


temperature. 


T 

c 

= 

critical temperature. 


T* 

= 

T/T . 
c 


AT* 

= 

(T-T )/T . 
c c 


T 

r 

= 

reference temperature (see eqs. (3.20) and (3.25). 


T 

= 

T/T . 
r 



Vi 



volume . 

AT*/|Ap*|^^^. 

-X at coexistence curve. 

(x+x )/x . 
o o 

power law exponent for specific heat (see Table I) . 

power law exponent for coexisting densities (see Table I) . 

power law amplitude for compressibility above T^ (see Table I). 

power law amplitude for compressibility below T^ (see Table I) . 

power law exponent for compressibility (see Table I) . 

power law exponent for chemical potential at T = T^ (see Table I). 

shear viscosity. 

lim ri. 
p o 

excess viscosity. 

normal viscosity. 

critical viscosity enhancement. 

parametric variable in linear model equation of state (see eq. (2.14)). 
correlation length. 

power law amplitude for correlation length above T^ (see Table I) . 
power law amplitude for correlation length below T (see Table I) . 

constant in eq. (3.11) for critical thermal conductivity 
enhancement . 

thermal conductivity. 

lim A. 

p-v o 

excess thermal conductivity. 

normal thermal conductivity. 

critical thermal conductivity enhancement. 

coefficients in eq. (3.17) for of carbon dioxide. 
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u 

y* 

Ay* 


y 


* 

c 


y*(T*) 

o 


y* 


V 


p 

Pc 

P* 


Ap* 


Xt 

x| 

y* 


= chemical potential. 

= yp^/p^. 

= {y(p,T) - yCp^,T)}p^/P^. 

= value of y* at critical point. 


= analytic background in the equation for the chemical potential 
in the critical region (see Table II) . 

= coefficients in eq. (2.19) for y*(T*). 

o 

= power law exponent for correlation length (see Table I). 

= density. 

= critical density. 

= p/p^. 

= (p-p^)/p^. 

= reference temperature (see eqs. (3.20) and (3.25)). 

= P/P^. 

= exponent in eq. (3.27) for critical viscosity enhancement. 

= p\- 

■ 


vili 



1 . 


Introduction 


Many thermophysical properties exhibit a singular behavior in the 
vicinity of the critical point. For instance, the isothermal compress- 
ibility, the thermal expansion coefficient and the specific heat of 
fluids all diverge at the critical point. Critical enhancement effects 
are also encountered in the behavior of the thermal conductivity and 
the viscosity near the critical point. It has always been difficult to 
formulate correlating equations that Incorporate accurately the singular 
behavior of the various thermophysical properties near the critical 
point. However, in recent years considerable progress has been made in 
the art of developing representative equations for the thermodynamic and 
transport properties of fluids in the critical region. The purpose of 
this technical report is to discuss these developments with an emphasis 
on practical applications. 

The critical point is a point of marginal thermodynamic stability. 

In the vicinity of the critical point, large-scale density fluctuations 
are present in the fluid. The spatial extent of these fluctuations can 
be characterized by a correlation length Near the critical point, 
this correlation length becomes much larger than the intermolecular 
interaction range. One can imagine that the critical enhancement effects 
are caused by the static and dynamic behavior of a system of clusters. 
Different states in the vicinity of the critical point correspond to 
different cluster sizes and the thermodynamic states can be interrelated 
by scaling laws. Moreover the properties of a system of clusters become 
independent of the nature of the individual molecules. This principle 
is commonly referred to as critical-point universality (Levelt Sengers 
et al. , 1977) . 

Recent developments have revealed that the simple scaling laws for 
the critical behavior of the thermodynamic properties become rigorously 
valid only in an asymptotically small region around the critical point, a 
region too small to be of interest in engineering applications. In order 
to formulate representations of the thermophysical properties that can 
be used in a larger range around the critical point one has two options. 

The first option, the more fundamental one, is to develop a theory for the 
higher order corrections to the asymptotic scaling laws. This approach 
is currently being developed (Levelt Sengers and Sengers, 1980). The 
other approach is to modify the asymptotic scaling laws using some semi- 
empirical adaptions. It is the latter approach which will be followed in 
this technical report. Hence, the emphasis is not on a test of the 
validity of the modern theory of critical phenomena, but rather on a 
presentation of some simplified procedures that enable one to characterize 
the behavior of the thermophysical properties of fluids in the critical 
region with reasonable accuracy. 

A comprehensive account of methods for describing equilibrium properties 
of fluids was presented in a previous NASA contractor report (Sengers and 
Levelt Sengers, 1977). The primary purpose of the present report is to ex- 
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tend the methods to a description of the transport properties of fluids 
in the critical region. However, the behavior of the transport properties 
near the critical point is closely related to the critical behavior of the 
thermodynamic properties and the correlation length. For this reason-, it 
was deemed desirable to include in this report a brief review of the 
behavior of the equilibrium properties in the critical region as well. 

There exists a close analogy between the behavior of the thermodynamic 
and transport properties of fluids near the gas-liquid critical point and 
those of binary liquids near the critical mixing point (Sengers and Levelt 
Sengers, 1977; Scott, 1978; Greer, 1978). This technical report deals only 
with the behavior of the properties of one-component fluids near the gas- 
liquid critical point. 

The part of this report concerned with the equation of state is based 
on research done in collaboration with the late T. A. Murphy. Representative 
equations for the thermophysical properties of carbon dioxide were formulated 
in collaboration with F. J. Cook and equations for the transport properties 
of steam in collaboration with J. T. R. Watson. 

Computer time for this project was provided by the Computer Science 
Center at the University of Maryland. 
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2 . 


Equilibrium properties 


2 . 1 Power laws 


Let A be the Helmholtz free energy, S the entropy, y the chemical 
potential per unit mass, the heat capacity at constant volume and Cp 
the heat capacity at constant pressure. We also find it convenient to 
introduce a symmetrized Isothermal compressibility defined as 


Xrj. 



( 2 . 1 ) 


The thermodynamic properties are made dimensionless by expressing them 
in terms of the critical temperature T^, the critical density p^, and the 
critical pressure . 

T* = T/T , p* = p/p , P* = P/P , A* = A/P V , y* = yp /P 
c c c c c c 


X*=X/^/P^ , S*=ST^/P^V 


C*=c T /P V , C*=C T /P V .(2.2) 
vvcc PPcc 


Note that the reduced extensive properties A*, S* and C* are all taken per 
unit volume, rather than per unit mass. The reason is ¥hat the singular 
part of the extensive thermodynamic properties per unit volume appear to 
be approximately S 3 anmetric or antisymmetric functions of p-Pj, (Vicentini- 
Missoni et al., 1969a; Levelt Sengers, 1974). As a consequence, Xt is also 
a symmetric function of p-p^, while the compressibility Rj, = p~i ( 8p/ 8P) .p is 
not (Sengers, 1973). We also define the reduced differences 

AT* = (T-T )/T , 

c c 

Ap* = (p-p )/p , (2.3) 

c c 

Ay* = [y(p,T) - y(p ,T)]p /P 
c c c 


To represent the singular thermodynamic behavior of fluids in the vicinity 
of the critical point, one first defines critical power laws. The ex- 
ponents of these power laws depend on the property considered and the path 
along which the critical point is approached. The special paths in the 
AT* - Ap* plane commonly considered are the coexistence curve or phase 
boundary (CXC) , the critical isochore (CIC) Ap* = 0 and the critical isotherm 
(CIT) AT*=0. These paths are indicated schematically in Fig. 1. The 
critical power laws defined along these paths are summarized in Table I. 

The critical exponents of the different thermodynamic properties are not 
independent. The scaling laws to be discussed below imply the exponent 
relations 


a = 2 - B(6+l) , Y = B(6-l) . (2.4) 
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Fig. 1 



Ap* 


CXC 


Special paths in the AT* versus Ap* plane for the definition 
of critical power laws. 
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It is generally believed that the critical exponents of fluids should 
be the same as those of spin systems represented by a Landau-Ginzburg 
Hamiltonian (Hubbard and Schofield, 1972; Nicoll, 1980). These critical 
exponent values have been calculated with considerable accuracy (Baker et al., 
1978; Le Guillou and Zinn-Justin, 1977, 1980). Experiments indicate 
agreement with these theoretical exponent values provided that the temperature 
of the gas-liquid critical point is approached to within a few 
hundredths of a percent (Levelt Sengers et al . , 1977). Efforts to 
develop a systematic theoretical description of the thermodynamic behavior 
over a larger range of temperatures and densities are currently in pro- 
gress (Ley-Koo and Green, 1977; De Pasquale et al., 1978; Nicoll and 
Chang, 1978; Levelt Sengers and Sengers, 1980). Here we continue to employ 
a more practical approach, adopted in the previous technical report, by 
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using the critical power laws over a larger range with effective exponent 
values deduced from the experimental data (Levelt Sengers and Sengers, 
1975; Levelt Sengers et al., 1976; Sengers and Levelt Sengers, 1977; 
Aharony and Ahlers, 1980). 


2.2 Scaling 


The critical power laws characterize the behavior of the thermo- 
dynamic properties along selected paths in the AT* - Ap* plane. We 
can generalize the description to the behavior of the thermodynamic 
properties at arbitrary temperatures and densities in the critical 
region by formulating scaling laws. For this purpose one defines a 
scaling variable x as 


X 


AT*/|Ap*|^'^^ 


(2.5) 


This scaling variable is chosen such that the coexistence boundary 
corresponds to a constant value of the scaling variable x: 


X = -X 

cxc o 



( 2 . 6 ) 


Instead of x, it is slightly more convenient to use the variable 


X + X 



Curves of constant y in the AT* - Ap* plane are indicated schematically 
in Fig. 2. The coexistence curve corresponds to y=0, the critical 
isotherm to y = 1 and the critical isochore to y=“. 
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AT* 



Fig. 2 Curves of constant y in the AT* versus Ap* plane. 


Sufficiently close to the critical point, the singular behavior 
of the thermodynamic properties, when reduced appropriately, becomes 
a function of the scaling variable y. To specify the thermodynamic 
properties in the critical region completely, the scaled contributions 
are supplemented with contributions that are analytic functions of the 
physical variables. The expressions thus obtained for a number of 
thermodynamic properties are listed in Table II (Sengers and Levelt 
Sengers, 1977). 
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Table II 


Scaled expressions for thermodynamic functions^ 


Chemical potential 

y*= y*(T*) + Ap*|Ap*|'^"^Dh(y) 

Compressibil ity 

= |Ap*|‘^~^D[6h(y) +^(l-y) ] 

Pressure 

P* = p* (T*) + D |^Ap*|Ap*|'^'^hCy) + |Ap*|'^^^|hCy) - a(y)}] 
Pressure coefficient 



Helmholtz free energy density 

A* = -P*(T*) + p*y*(T*) + |Ap*|^'^^Da(y) 


Entropy density 


dP’' 

o 

“df* 


- 


dy*(T*) 

o 

dT* 


dy 


Heat capacity 

C* d2p* d^y*(T*) 


dT* 


dT* 


-a/B n . 

|4p,| JLi^ 

X dy 


P*(T*) and y*(T*) are analytic functions of T* . 
o o 
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Here the singular behavior of the thermodynamic properties is expressed 
in terms of two scaling functions h(y) and a(y) which are related by the 
differential equation (Griffiths, 1967; Levelt Sengers et al., 1976): 

Bh(y) = B(6+l)a(y) + (1 - y)^|^ . (2.8) 

The function h(y) is normalized such that h(l) = 1. 

As an example, we consider the scaling law for the chemical potential 


Ay* 


Ap* I Ap* 


6-1 


Dh(y) 


(2.9) 


Normally, the chemical potential is a function of two variables, density 
and temperature. However, the scaling law (2.9) implies that when "scaled" 
chemical potential data Ay*/DAp* | Ap* | <S-1 are plotted as a function of the 
"scaled" temperature y, the data should collapse onto a single curve. As 
an example, we show in Fig. 3 such a plot for steam. This plot is based 
on the experimental data of Rlvkin and coworkers (1962, 1963, 1966) as 
analyzed by Levelt Sengers et al . (1976). 

2 . 3 Equations of state for the critical region 

2.3.1 MLSG equation of state 

The equation of state is commonly defined as a functional 
relationship between the pressure P, the volume V and the temperature T. 
However, in the critical region the density p is a more suitable variable 
than the volume V (Vicentini-Mlssoni et al., 1969a). Since y is the 
thermodynamic function conjugate to the density p, the equation of state 
considered here is the functional relationship between the chemical 
potential y, the density p and the temperature T. Thus to specify the 
scaled equation of state, we need to specify the function h(y) in (2.9). 

One equation, proposed for this purpose, is the so-called NBS- or MLSG- 
equation (Vicentini-Mlssoni et al., 1969a, b; Levelt Sengers et al., 1976; 
Sengers and Levelt Sengers, 1977) 


h(y) 


y 



(Y-l)/2g 


( 2 . 10 ) 
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Xo 


Fig. 3 Scaled chemical potential data for steam as a function of 

y= (x + Xq)/xq. Reduction parameters are taken from Table IV. 


where E is a constant. As a consequence, the symmetrized isothermal 
compressibility X* is represented by 


X 


* 

T 


Ap* 


-y/3' 


[(1 + E)/(1+Ev^^)1 
D[6+ (y-l)(6- 3"^ 


(Y-D/26 

+ Ey ^^)/(1 + Ey ^^)] 


( 2 . 11 ) 


Note that at the critical isochore Ap*=0 and eq. (2.11) reduces to 


= r|AT*r^ (2.12) 

In terms of the MLSG equation of state, the amplitudes B, F and F' of the 
critical power laws, defined in Table I for the coexistence curve and the 
compressibility, are related to D and x^ by 
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B = 

0 

(2.13a) 


(2.13b) 

r' = -^(i + eI 

(2.13c) 

2.3.2 Linear model equation of state 


The MLSG equation has the advantage that it enables one to 
calculate the compressibility directly for given values of AT* and 

Ap*. However, it has the disadvantage that the expression (2.10) for 
h(y) cannot be integrated analytically to yield an explicit expression 
for a(y) and, hence, for the pressure. This problem is solved by 
introducing two auxiliary ("parametric") variables r and 0 instead of 
the variables AT* and Ap* . 

Here we consider the linear model, originally 
Schofield (1969). It is defined by 

introduced by 


AT* = r(l-b^e^) 

(2.14a) 

Ap* = kr^0 , 

(2.14b) 

Ay* = ar^'^ed - 0^) 

(2.14c) 


where a, k and b are constants. This model is chosen such that 0= ±1 
correspond to the two branches of the coexistence curve, 6 = ± 1/b to the 
critical Isotherm and 0 = 0 to the critical isochore. The resulting 
parametric expressions for a number of thermodynamic properties are given 
in Table III. The amplitudes of the critical power laws defined in Table I 
are related to the linear model constants by (Hohenberg and Barmatz, 1972; 
Sengers and Levelt Sengers, 1977) 


k/(b^-l)^ 

(2.15a) 

a(b^-l)b'^"^/k^ 

(2.15b) 

k/a , 

(2.15c) 

(b^- 1) { l-b“ (1-2B) f k/2a 

(2.15d) 
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(?..15e) 


a"*” = ak(2-a) (l-a)a 
A~ = ak(2-a)(l-a)a(p^+ P^+ P^)/b^-l)^""“ . (2.l5f) 


Table III Linear model equations for thermodynamic properties 


Variables 


AT* = r(l-b^e^D 


Ap* = r^k6 


Chemical potential 


y* = y*(T*) + r^\e(l - 0^) 
o 


Compressibility 


X* = ^ 

^T a 


l-(l-2g)b^6^ 


1-30^ + 1 + 266ci-e^)| 


Pressure 


p* = p*fT*) + r^^ae(l-e^) + r^*^^‘^^^akp(6) 


Pressure coefficient 


fip*C7’*'l 2 2 

1-a , ^ 36-1 on 6fl-0 )-a-30 ) 

— + r aks(6) + p*r^ agO — ' ■ y-y 

* l-(l-2g)b 0 


/9P* \ ^ ^ 
V3T*/ d' 


Helmholtz free energy density 


A* = - p*(T*) + p*y*(T*) + “ak{0^(l-0^) - p(0)} 

o o 


Entropy density 


S* = 


dP*(T*) 

dT* 


y*(T*) 

- p* “aks(0) 
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Heat capacities 


* 


C 


V 


T* 


d^P*fT*) 

2 

dT*"^ 


d^y*(T*) 

2 

dT* 


^ -a 1 

+ r ak 


II 

1 



g)sCe)-3es'fQl 

(l-23)b^e^ 






A computer program calculating the parametric variables r and 0 for given 
temperatures and densities is presented in the Appendix. A similar program 
has been published by Holdover (1978) . 

2 

The linear model constant b is sometimes identified with (Schofield 
et al . , 1969) 


h2 = <S-3 

(6-1) (1-26) 


(2.16) 


The linear model with condition (2.16) is referred to as the restricted 
linear model. Further details can be found in the previous technical report 
(Sengers and Levelt Sengers, 1977). 

2 .4 Universality 

The principle of universality asserts that the critical exponents and the 
scaling function h(y) should be the same for all 3-dimensional systems with 
short range forces. The only system dependent constants in the scaling laws 
are the coefficients D and Xq which specify the amplitude of the critical 
isotherm and the coexistence curve, respectively (Levelt Sengers, 1977a). In 
terms of the MLSG equation of state and the linear model equation of state, 
this principle requires the constants E and b^ to be universal. 

Since we are applying the scaling laws outside their asymptotic range using 
effective critical exponent values and effective values for E and b^ , the 
principle of universality will only apply approximately. In order to Illustrate 
to what extent this empirical extension of the universality principle is 
adequate, we reproduce in Fig. 4 scaled chemical potential data for ^He, ^He, 

Xe, CO 2 and H 9 O) (Levelt Sengers, 1974; Levelt Sengers et al., 1976; Sengers and 
Levelt Sengers, 1977). The data can be represented approximately by a single 
function of the scaling variable y= (x+Xq)/xq. The data cover a range in 
temperatures and densities bounded by 

5 X 10 "‘^S|AT*| <3 X 10 "^ , ]Ap*|s0.25 . (2.17) 


The principle of corresponding states would imply identical behavior when 
the thermodynamic properties are reduced using the critical parameters T^, , , 

Pq. The fluids included in Fig. 4 do not satisfy corresponding states. 
However, the universality hypothesis is a generalization of corresponding 
states, saying that identical behavior of the singular critical behavior of 
the thermodynamic properties is obtained by use of two additional constants, 

D and x^, or in terms of the linear model, a and k. In Table IV, taken from 
the previous technical report (Sengers and Levelt Sengers, 1977), we have 
listed these parameters for a number of fluids. In this table we have added 
parametric values for isobutane, taken from an NBS Internal Report (Waxman 
et al., 1978). 
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4 
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1 I I 

10 10 ^ 10 ^ 
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Fig. 4 Scaled chemical potential data for ^He, ^He, Xe, CO 2 and H 2 O as 

a function of y = (x + Xq)/xq. Reduction parameters are taken from 
Table IV (Sengers and Levelt Sengers, 1977). 


15 



Table IV. 

Critical region parameters for 

1 number of fluids with 



effective 

universal 

exponent values (Sengers and Levelt 

Sengers . 

1977) . 






P 

c 

MPa 

Pc 

kg/m^ 

T 

c 

K 

X 

0 

D 

k 

a 

c xio^o 

0 

m 

^He 

0,11678 

41.45 

3.3099 

0.489 

3.16 

0.924 

4.58 

2.7 


0.22742 

69.6 

5.1895 

0.369 

2.86 

1.021 

6.40 

2.2 

Ar 

4.865 

535. 

150.725 

0.183 

2.43 

1.309 

16.1 

1.6 

Kr 

5.4931 

908. 

209.286 

0.183 

2 .43 

1.309 

16.1 

1.7 

Xe 

5.8400 

1110. 

289.734 

0.183 

2.43 

1.309 

16.1 

1.9 

p-H^ 

1.285 

31.39 

32.935 

0.260 

2.50 

1.156 

9.6 

1.9 

N 

2 

3.398 

313.9 

126.20 

0.164 

2.32 

1.361 

18.2 

1.6 

°2 

5.043 

436.2 

154.580 

0.183 

2.36 

1.309 

15.6 

1.6 

«2° 

22.06 

322.2 

647.13 

0.100 

1.28 

1.622 

21.6 

1.3 

D 2 O 

21.66 

357. 

643.89 

0.100 

1.28 

1.622 

21.6 

1.3 

CM 

0 

7.3753 

467.8 

304.127 

0.141 

2.16 

1.436 

21.3 

1.6 

NH 3 

11.303 

235. 

405.4 

0.109 

1.47 

1.573 

21.4 

1.4 


3.7605 

730. 

318.687 

0.172 

3.06 

1.337 

23.9 

2.0 

CH. 

4 

4.595 

162.7 

190.555 

0.164 

2.17 

1.361 

17.0 

1.7 

^2»4 

5.0390 

215. 

282.344 

0.168 

2.32 

1.350 

17.5 

1.9 

S»6 

4.8718 

206.5 

305.33 

0.147 

2.17 

1.416 

20.2 

1.8 

S»8 

4.247 

221. 

369.82 

0.137 

1.96 

1 .451 

20.2 

2.0 


3.631 

227. 

407.85 

0.140 

2.07 

1.441 

20.8 

2.2 

Notes ; 

a = 0. 100, 

B= 0.355, Y= 1-190, 6 

= 4.352, V = 

0.633, 

E = 0.287 

b^ = 1.3909 
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2 . 5 A fundamental equation for the critical region 


The scaled equations introduced in the previous section enable us to 
calculate the singular critical behavior of the thermodynamic properties. 
However , in order to obtain the complete thermodynamic behavior , we need to 
specify the analytic background functions P*(T*) and y*(T*) (see Tables II 
and III). This is accomplished by expanding P*(T*) and yg(T*) in a Taylor 
series in terms of AT* (Levelt Sengers, 1977b). 


P*(T*) = 1 + Z P*(AT*)^ , (2.18) 

° i=l ^ 

y*(T*) = U* + z: y*(AT*)^ . (2.19) 

o c . - 1 

1=1 

Based on linear model parameters earlier determined by Murphy et al., 
(1973, 1975), these background parameters have been determined for steam and 
carbon dioxide. For steam this task was done by Levelt Sengers (1977b) and 
the fundamental equation parameters are reproduced in Table V. For carbon 
dioxide this task was accomplished by Cook and Sengers as part of this project. 
The parameters for CO 2 , deduced from the experimental data of Michels et al. 
(1937, 1948, 1952), including a correction for the temperature scale (Levelt 
Sengers and Chen, 1972; Levelt Sengers et al., 1976), are also given in Table 
V.t 


2 . 6 Correlation length 

We shall try to relate the critical enhancement of the transport 
properties to the thermodynamic properties and the correlation length 5 
which measures the sixe of the critical fluctuations. For this purpose we 
calculate the correlation length using an approximate relationship with the 
compressibility introduced in the previous technical report (Sengers and 
Levelt Sengers, 1977) 


e , (2.20) 

where the correlation function exponent v is related to the thermodynamic 
critical exponents by the hyperscaling relation 


In this report all 
Practical Scale of 


temperatures are quoted in terms of the International 
1968, referred to as IPTS '68 (Douglas, 1969). 
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Table V Fundamental equation parameters based on 
the linear model with effective exponent values 


steam 


t 


carbon dioxide 


ft 


Critical parameters 


p = 324.383 kg/m 
c 

P = 22.0477 MPa 
c 

= 647.073 K (IPTS '68) 


= 467.8 kg/m 

P = 7.37516 MPa 
c 

= 304. 127 K ( IPTS '68) 


Critical exponents 


a = 0.08712 
e = 0.3505 

Y = 1.2119 
6 = 4.4576 

V = 0.6367 


a = 0.128 
B = 0.3486 

Y = 1.175 
S = 4.37 

V = 0.624 


Linear model parameters 


a = 24.0999 


a = 24.4878 


k = 1.6837 


k = 1.745 


b = 1.4099 


b = 1.70 


Coefficients F* 
1 


P^ = 7.87425 

P* = -25.8448 
2 

Pt = 3.96522 


6.9992 


= -8.8535 


Coefficients uf 
1 


Correlation length 
amplitude 


Range 


p* = -34.05 
p* = -16.59 
p* = 59.84 


1.31 


^2 

U* 


13.0231 

-49.2265 

423.428 


e = 1.5 X 10“^°m 
o 


[at* I S 0.05, I Ap* I <0.25 | AT* | < 0.03, |Ap*| S 0.28 


From Levelt Sengers (1977b) . 

From Cook and Sengers (unpublished report, 1974). 
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V = B(6+ l)/3 . (2.21) 

The amplitude for the correlation length can he related to the critical power 
law amplitudes B and F, defined in Table I, by 

5 = R(Fk T /B^P (2.22) 

O B C C 

where R is a universal constant whose value is approximately equal to 0.7 

(Sengers and Levelt Sengers, 1977; Sengers and Holdover, 1978). Values 

for the correlation length amplitude ^ are included in Tables IV and V. 

o 

The above equations allow us to calculate the correlation length ^ from 
the symmetrized reduced compressibility X'* by the use of either the MLSG 
equation (2.11) or the linear model equation for xf given in Table III. 

3 . Transport properties 

3 . 1 Introductory remarks 

The thermal conductivity X and the shear or dynamic viscosity iq exhibit a 
critical enhancement in the vicinity of the critical point. For a survey of 
the history of the subject the reader is referred to some earlier reviews 
(Sengers, 1966, 1971, 1972, 1973). 

In Fig. 5 we show schematically the thermal conductivity X as a function 
of the density p at three temperatures Tq^ < T£ < T^ above the critical temperature 
Tj, . In order to describe the critical enhancement, we separate X into a normal 
or background thermal conductivity X and a critical thermal conductivity 
enhancement AX (Sengers and Keyes, 1971). 

Outside the critical region, the_thermal conductivity X is to be_identlf led 
with the normal thermal conductivity X. Inside the critical region, X is defined 
empirically by extrapolating the behavior of the normal thermal conductivity 
outside the critical region. The normal thermal conductivity X may in turn be 
written as 


X(p,T) = X (T) + X (p,T) , (3.2) 

o e 

where Xq(T) is the thermal conductivity in the limit of low densities at the same 
temperature and where Xg(p,T) is commonly referred to in the engineering 
literature as the excess thermal conductivity. The excess thermal conductivity 
is a slowly varying function of temperature and increases monotonically with 
density (Diller et al . , 1970; Hanley et al., 1972; Sengers, 1968). 
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Fig. 5 


Schematic representation of the thermal conductivity A(p,T) as a 
function of the density p at three supercritical temperatures 


The behavior of the viscosity p as a function of density at a temperature 
slightly above the critical temperature is shown schematically in Fig, 6, In 
analogy to (3.1) and (3.2) we separate the viscosity p into a normal viscosity 
n and a critical viscosity enhancement An 


n = n + An 


(3.3) 


with 


n(p,T) = n^(T) + ng(p,T) . (3.4) 

However, there are some qualitative differences between the critical enhancements 
of the thermal conductivity and of the viscosity. First, the thermal con- 
ductivity enhancement is a "strong" enhancement extending over a large range of 
temperatures and densities, while that of the viscosity is a "weak" enhancement 
limited to a small range of temperatures, typically a few degrees, around the 
critical temperature (Sengers, 1973). Secondly, theoretical and experimental 
evidence indicates that the critical viscosity enhancement is a multiplicative 
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Fig. 6 Schematic representation of the viscosity r|(p,T) as a function of 

the density p at a temperature slightly above the critical temperature. 


enhancement, l.e. At) is itself proportional to the normal viscosity (Kawasaki, 

1971, 1976; Perl and Ferrell, 1972; Sengers 1971, 1973; Ohta, 1977; Calmettes 
1977). Thus, in the case of the viscosity, we prefer to consider the relative 
critical enhancement Ar)/r| or the viscosity ratio n/r). 

Another transport property of interest is the thermal diffusivity 
D.J. = A/pCp, where Cp is the specific heat at constant pressure per unit mass. It 
can be determined experimentally from measurements of the decay rate of the 
critical fluctuations with the use of light scattering techniques (Mountain, 

1966; Benedek, 1968; McIntyre and Sengers, 1968; Cummins and Swlnney, 1970; 

Fleury and Boon, 1973) . This method has the advantage that the thermal diffusivity 
is measured while the fluid remains in thermodynamic equilibrium but it is in 
practice limited to a small range around the critical point, where the scattered 
light intensity is sufficiently large for the method to be applicable. It is 
also possible to determine the thermal diffusivity from the temperature field 
in the fluid near a heated boundary, made visible by holographic interferometry 
(Becker and Grigull, 1975, 1978; Schmitt and Beer, 1977). Although the thermal 
conductivity diverges at the critical point, the specific heat at constant 
pressure Cp increases much faster, so that the thermal diffusivity goes to zero 
at the critical point. 
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As mentioned in the previous technical report (Sengers and Levelt Sengers, 
1977) , the modern theory of critical phenomena is based on a renormalization- 
group approach. For a more comprehensive introduction to the subject the 
reader is referred to the literature (Barber, 1977; Bomb and Green, 1976; 

Pfeuty and Toulouse, 1977; Wallace and Zia, 1978; Wilson and Kogut, 1974). 
Subsequently, the method has been extended to dynamic critical phenomena as 
well (Siggia et al., 1976; Kawasaki and Gunton, 1976; Gunton and Kawasaki, 

1976; Hohenberg and Halperin, 1977; Ernst, 1978; De Dominici and Peliti, 1978; 
Mazenko, 1978). Another approach to the theory of dynamic critical phenomena 
is commonly referred to as mode-coupling theory. In this approach the time 
correlation function expressions for the transport coefficients are evaluated 
by use of the methods of generalized hydrodynamics and the critical anomalies 
arise from non-linear coupling between hydrodynamic modes (Kadanoff and Swift, 
1968; Kawasaki, 1970, 1971, 1976; Pomeau and Reslbois, 1975). For a discussion 
of the relationship between the mode-coupling theory and the dynamic renormal- 
ization-group theory the reader is referred to a recent survey presented by 
Gunton (1979). These theories yield expressions for the transport coefficients 
near the critical point as a function of the wavelength of the fluctuations. In 
this technical report we only consider the hydrodynamic transport coefficients 
corresponding to the long wavelength limit. 


The dynamic renormalization-group theory predicts that the thermal con- 
ductivity and the viscosity near the critical point should diverge as (Siggia 
et al., 1976; Hohenberg and Halperin, 1977) 


k T 

AA = A Pc 

birriC p 


(3.5) 


An = n “ » (3.6) 

where kg is Boltzmann's constant and where A is a universal constant and (}) a 
universal exponent independent of the nature of the fluid. 


The mode-coupling theory leads to sets of coupled integral equations which 
are difficult to evaluate. In first approximation one finds (Kawasaki, 1970, 
1971, 1976; Ferrell, 1970; Perl and Ferrell, 1972) 

kpl k^T 

~ 6Trri5 ^^p ~ 6irriC ^^*^p ’ 

^ Jln(qO 

^ 15F 


(3.7) 

(3.8) 
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where q is a slowly varying function of temperature (Oxtoby and Gelbart, 1974), 
but which will be treated here as an adjustable constant. Equation (3.7) may 
be interpreted physically by assuming that the critical part AX/pc^ of the 
thermal diffusivity D-j is determined by the Stokes-Einstein diffusion of clusters 
with radius 5 • 


(3.6) 


Since Ar|/r|«l, equation (3.8) may be written in the form of a power law 


= 1 + Jin(qO‘^ (qO 


(3.9) 


with (j) == 8/15 tt^ = 0.054 (Ohta, 1975; Ohta and Kawasaki, 1976). The exponent (J) has 
also been calculated from a perturbation expansion up to second order in 
e=4-d, where d is the dimensionality of the system, yielding <})=: 0.065 (Siggia 
et al., 1976; Gunton, 1979). Expanding A up to first order in e, Siggia et al. 
estimated A- 1.2. On comparing (3.7) with (3.5) we note that the mode-coupling 
theory implies in first approximation the value A= 1.0. Recent experimental 
studies very close to the critical point have indicated that the value of A is 
Indeed close to unity (Burstyn et al., 1980; Giittinger and Cannell, 1980). 


As mentioned in Section 2, the asymptotic scaling laws for the equilibrium 
properties become valid rigorously at temperatures within a few hundredths of a 
percent from the critical temperature. Due to the additional approximations in 
the derivation of the asymptotic power laws for the transport properties, the 
range of validity of equations Introduced above will be even smaller. Just as 
in the case of the thermodynamic properties, we adopt here a practical approach 
in which the thermal conductivity and the viscosity in the critical region are 
presented by (3.5) and (3.9), but allowing for effective values of the amplitude 
A and the exponent (}) which may differ slightly from the values predicted 
theoretically. 


3 . 2 Thermal conductivity 

3.2.1 Thermal conductivity equation 

The asjrmptotlc equation (3.5) for the critical thermal conductivity 
enhancement will only be applicable in a relatively small region of temperatures 
and densities around the critical point. However, the actual observed critical 
thermal conductivity enhancement covers a large range of temperatures and 
densities. For instance, in the case of carbon dioxide a critical enhancement 
is noted in a temperature range of 70°C above the critical temperature 
(Le Neindre et al., 1973). In practice, one needs a representative equation 
covering the entire range of the critical enhancement. For this purpose we 
adopt an equation of the form (Hanley et al., 1974b, 1976). 


X = X + AX 


(3.10) 


with 
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(3.11) 


k T 

AA = A P(c -c )F(AT*,Ap*) 

6pn5 p V 

Here F(AT*,Ap*) is an empirical crossover function such that 

lim F(AT*,Ap*) = 1 , lim F(AT*,Ap*) = 0 . (3.12) 

|AT* 1^0, |Ap* |-^ 0 I AT* !■>«>, I Ap* ® 

In practice, we continue to use a crossover function earlier introduced by 
Hanley et al. (1974b, 1976) : 

F(AT*,Ap*) = p*" exp[-{A^(AT*)^ + B^(Ap*)^}] . (3.13) 

The crossover function insures that equation (3.11) for AA reduces to the 
theoretically predicted equation (3.5) upon approaching the critical point, 
while A reduces to the normal thermal conductivity A far away from the critical 
point . 

Equation (3.11) relates the critical thermal conductivity enhancement AA 
to the thermodynamic properties, the viscosity r| and the correlation length 
Using the thermodynamic relationship between the specific heat at constant 
pressure and constant volume 


p(c - c ) = 
P V 


P 

c 

T 

c 




(3.14) 


and the relationship (2.20) between the correlation length and the symmetrized 
compressibility xf > we may rewrite (3.11) as 


with 




(3.15) 




v/y 


(3.16) 


Since the critical viscosity enhancement is much weaker than the critical thermal 
conductivity enhancement, we may approximate the viscosity n in (3.11) and (3.15) 
by the normal viscosity r|. This approximation is adequate within the accuracy 
with which the thermal conductivity near the critical point is known experi- 
mentally. 
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The function F In (3.11) and (3.15) is an empirical crossover function and 
the parameters and as well as the exponent n, in (3.13) are to be 
determined from the experimental thermal conductivity data. As a first approx- 
imation, these parameters were treated by Hanley et al. (1974b, 1976) as con- 
stants independent of the nature of the fluid. This assumption is strictly 
only valid for noble gases, where the thermal conductivity can be expected to 
satisfy corresponding states. In practice, these parameters, as well as the 
optimal form of the crossover function, may vary somewhat from fluid to fluid. 

The calculation of the thermal conductivity enhancement from (3.15) requires 
knowledge of the symmetrized isothermal compressibility x* and of the pressure 
coefficient (3P*/3T*)p. Since the compressibility is strongly divergent, it is 
imperative that xf calculated from an appropriate scaled equation. For a 
large number of fluids the compressibility xf can be calculated from the para- 
meters quoted in Table IV. However, in order to calculate (3P*/3T*)p an 
equation for y*(p*,T*) is insufficient and a more complete equation of state is 
required. A more complete equation of state based on the scaling laws is 
currently available for few substances only, namely for steam (Levelt Sengers, 
1977b) and for carbon dioxide (Basu and Sengers, 1980) as discussed in Section 
2.5, and also for Isobutane (Waxman et al., 1978). However, since (3P*/3T*)p 
remains finite at the critical point, already a substantial improvement is 
obtained when Xt calculated from the parameters presented in this report, 
even while (3P*/3T*)p is calculated from a conventional equation of state 
(Hanley et al., 1976). 

3.2.2 Thermal conductivity of carbon dioxide 

In order to test the applicability of the proposed thermal conductivity 
equation, we consider a few representative examples where adequate experimental 
information is available. In particular for carbon dioxide, the thermal con- 
ductivity has been determined over a wide range of temperatures and densities 
around the critical point. The thermal conductivity in the critical region up 
to 75°C was measured by one of the authors in collaboration with Michels and 
van der Gullk using a parallel plate method (Michels et al., 1962). The 
temperature range was subsequently extended by Le Nelndre et al. (1973) up to 
700°C using a concentric cylinder method. The thermal diffusivity in the 
critical region of CO 2 was measured accurately by Becker and Grigull (1978) 
from a holographic interferometric analysis of the transient temperature field 
in the gas under a heated horizontal plate. 

The data_of Le Nelndre et al. (1973) enable us to present the normal thermal 
conductivity X as 

, (3.17) 


A = 


A (T) 
o 


3 

+ E 
1=1 


with 
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(3.18) 


X^(T) = Z a,^/T’" 

° k=0 ^ 


Here T and p are dimensionless variables defined as 


T = T/T^ . p = p/p^ , (3.19) 

with respect to a reference temperature Tj. and a reference density pj.. In 
the case of CO 2 we adopt, for convenience only, the reference values 

= 100 K , p^ = 100 kg/m^ . (3.20) 


The coefficients aj^^ in (3.18) were determined by Watson (1976b) and are 
reproduced in Table VIA. The coefficients in (3.17) were presented in an 

earlier technical report (Sengers, 1972) and are reproduced in Table VIB. 


Table VI. 


Coefficients in the equation for the normal thermal conductivity 
X of C0„ 


VIA The coefficients a. 


X 




32.3537 

m.K/W 


77.1034 

m.K/W 


618.115 

ra.K/W 

4 = - 

613.848 

m.K/W 


VIB The coefficients X. 

1 




-9 

-1 ,-i 

= 0.33789 

X 

10 

Wm K 



-3 

-1 -1 

X 2 = 0.27725 

X 

10 

Wm K 



-4 

-1 -1 

X^ = 0.47396 

X 

10 

Wm K 


Using the parameter values 


n = 0 , A = 1.02 , A^ = 39.8 


Bx = 5.45 


(3.21) 
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we present in Fig. 7 a comparison of the proposed thermal conductivity 
equation with the thermal conductivity data of Michels et al. (1962) and 
those reported by Becker and Grlgull (1978) . For this purpose we adopted 
the values n= 0 and A= 1.02, while and were treated as adjustable 
parameters (Basu and Sengers, 1980). The thermodynamic properties in a 
region bounded by [AT*] S0.03 and |Ap*j S0.28, and also the constant 5 , 
were calculated using the parameter values given in Table V. Outside ° 
this region, we used tabulated thermodynamic property data values based on 
the experimental data of Michels and coworkers (Michels and Michels, 1935, 
1937; Michels et al., 1937; Michels and de Groot, 1948; Michels and Strljland, 
1952) . The viscosity was calculated from an equation formulated by Watson 
(1976a) . 



Fig. 7 Thermal conductivity data for carbon dioxide in the critical region 
along isotherms as a function of density. The curves indicate the 
values calculated from the representative equation proposed in this 
report . 
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It is noted that the property determined experimentally by Becker and 
Grigull (1978) is the thermal diffusivity and not the thermal conductivity. 

Hence, a comparison between the two data sets is affected by the accuracy 
with which the specific heat c close to the critical point can be calculated. 

The thermal conductivity data, ^reported by Becker and Grigull and reproduced 
in Fig. 7, were deduced from the experimental thermal diffusivity data, when 
the specific heat was calculated from the fundamental equation presented in 
Section 2.5 with 1^= 304.127 K as determined by Holdover (1974). However, 
the actual value of observed by Becker and Grigull in their experiment was 
304.144 K. If the latter value of is used, the thermal conductivities 
deduced from the thermal diffusivity data become slightly larger close to the 
critical point, leading to an almost perfect agreement with the data of 
Michels et al. (1962) at temperatures down to 305.21 K, but yielding a somewhat 
larger peak at 304.32 K. In view of the fact that the peak heights at 304.32 K 
and 305.21 K are sensitive to the value adopted for Tj,, we conclude from Fig. 7 
that our equation yields a reasonable representation of the thermal conductivity 
In the critical region. 

3.2.3 Thermal conductivity of steam 

The existence of a critical enhancement in the thermal conductivity 
of steam was demonstrated experimentally by Le Neindre et al. (1973) with a 
concentric cylinder apparatus and by Sirota, Latunin and Belyaeva (1973) with 
a parallel plate apparatus. Subsequently, Sirota et al. (1974, 1976) have 
published a large amount of data for the thermal conductivity of steam in the 
critical region. 

In order to represent the experimental thermal conductivity data in the 
critical region by the equations proposed in Section 3.2.1, we need again an 
equation for the normal thermal conductivity A = A^(T) + Xg(p,T) . In the case of 
carbon dioxide discussed in the previous section, this task was considerably 
simplified by the observation that the excess thermal conductivity Ag(p,T) is 
independent of the temperature T in a large range of temperatures outside the 
critical region (Le Neindre et al., 1973). However, the assumption that the 
excess thermal conductivity Ag(p,T) is independent of the temperature over a 
large range of temperatures, is not justified in the case of steam. For this 
reason we do not want to use equation (3.17) with values for the coefficients 
Aj^ deduced from thermal conductivity data far away from the critical temperature. 
In our first approach, we retained equation (3.17) but determined the coefficients 
A^ simultaneously with the parameters in the equation (3.15) for AA from the 
experimental data in the critical region. This procedure yields a satisfactory 
representation of the experimental thermal conductivity data of steam in the 
critical region (Basu and Sengers, 1977). However, the equation thus obtained 
cannot be readily extrapolated to represent the thermal conductivity of water 
and steam over a large range of temperatures and densities. 

In order to accomplish the latter purpose, we adopted for A an equation of 
the form 

- ^ 5 , . _ . 

A = A (T)exp[p S Z bC. (i- 1) (p - 1)^ ] , (3.22) 

° i=0 3=0 f 
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with 


\ (T) = E 

° k=0 


(3.23) 


Here T and p are again dimensionless variables defined as (cf. (3.19)) • 


T = T/T^ , p = p/p^ . (3.24) 

For steam we took the reference values 

= 647.27 K , p^ = 317.763 kg/m^ . (3.25) 

The choice of this equation for X was motivated by the observation that an 
equation of the same form yields an accurate representation of the viscosity 
of water and steam in a large range of temperatures and densities outside the 
critical region (Hendricks et al., 1977; Nagashima, 1977; Watson et al. , 1980a, b) . 
The reference values, defined in (3.25), are those currently used in the 
international formulations adopted by the International Association for the 
Properties of Steam (Nagashima, 1977; Schmidt, 1979; Meyer et al., 1980). These 
reference values are close to, but not equal to, the critical temperature and 
density of steam (Balfour et al. , 1980). The parameters a^ in equation (3.23) 
for Xq were earlier determined by Aleksandrov and Matveev (1976) ; they are 
reproduced in Table VIIA. A comparison of the equation for Xq with the 
experimental thermal conductivity data of water vapor and steam at low densities 
was presented in a previous technical report (Sengers and Basu, 1977). The 
coefficients b^^ In equation (3.22) were determined by Watson (1978) in 
collaboration with the present authors; they are listed in Table VIB. Details 
of how these coefficients were obtained, as well as a comprehensive comparison 
of the equation with all available experimental data is presented elsewhere 
(Sengers et al., 1980). 

To represent the thermal conductivity of steam in the critical region, we 
combine (3.22) with (3.13) and (3.15) using the parameter values 

n = 1/2 , A = 1.20 , A^ = 18.66 , B^ = 1.00 , (3.26) 

as adopted in an earlier analysis (Basu and Sengers, 1977) and using the critical 
region parameters presented in Section 2.5. The parameters for the fundamental 
equation for the thermodynamic properties of steam in the critical region, listed 
in Table V, enable us to calculate these properties in a region' bounded by 
|AT*| so. 05 and |Ap*fs0.25. Outside this region we used the 1968 IFC Formulation 
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Table VII Coefficients in the equation for the 
normal thermal conductivity X of steam 


VIIA The coefficients aj^ 


a^ = +2.022 23 m.K/W 
0 

a^ = +14.111 66 m.K/W 
= +5.255 97 m.K/W 
a^ = -2.018 70 m.K/W 


VIIB 


The coefficients b. . 

ij 


a 

0 







2 



3 




4 

1 

1 

0 

+1.329 

304 

6 

+1. 

701 

836 

3 

+5 

224 

615 

8 

+8.712 

767 

5 

-1. 

852 599 9 


-0.404 

524 

37 

-2. 

215 

684 

5 

-10 

124 

111 


-9.500 

061 

1 

+ 0. 

934 046 90 


+0.244 

094 

90 

+ 1. 

651 

105 

7 

+4 

987 

468 

7 

+4.378 

660 

6 

0.. 

0 


+0.018 

660 

751 

-0. 

767 

360 

02 

-0. 

272 

976 

94 

-0.917 

837 

82 

0. 

0 

4 

-0.129 

610 

68 

+0. 

372 

833 

44 

-0 

430 

833 

93 

0.0 



0. 

0 

5 

+0.044 

809 

953 

-0. 

112 

031 

60 

+0. 

133 

338 

49 

0.0 



0. 

0 

J 
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for Scientific and General use (IFC, 1968; Rivkin et al., 1978). The viscosity 
was calculated from the international formulation adopted by the International 
Association for the Properties of Steam (Hendricks et al., 1977; Nagashlma, 

1977) . The experimental thermal conductivity data were originally obtained by 
Sirota et al. (1974, 1976) as a function of pressure. The densities were then 
calculated from the experimental pressures and temperatures with the use of the 
equation of state data. However, the maxima in the thermal conductivity isobars 
as a function of density as reported by Sirota et al. do not seem to approach 
the critical density as observed for other fluids such as CO 2 . The analysis 
does depend on the value adopted for Tj, which was not determined in the experi- 
ments of Sirota et al. We, therefore, recalculated the densities, treating T^ 
as an adjustable parameter and selecting 1^= 646.923 K as discussed elsewhere 
(Basu and Sengers, 1977). A comparison of our representative equation with the 
experimental thermal conductivity data at the recalculated densities is shown 
in Fig. 8. 



DENSITY, kq/m’ 


Fig. 8 Thermal conductivity data for steam in the critical region at four 
different isobars, measured by Sirota et al. (1974, 1976). The 
curves indicate the values calculated from the representative equation 
proposed in this report. 
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The equation for the thermal conductivity of steam presented here was 
adopted by the International Association for the Properties of Steam as the 
recommended interpolating equation for scientific use (Kestin, 1978, 1980; 
Scheffler et al. , 1979; Schmidt, 1979; Meyer et al., 1980). Further details 
will be presented in a forthcoming publication (Sengers et al., 1980). 

3.3 Viscosity 


3.3.1 Viscosity equation 


The asjraiptotic power law (3.9) for the critical viscosity enhancement 
is expected to become valid at temperatures and densities sufficiently close to 
the critical point such that qC » 1. In order to describe the critical 
viscosity enhancement more completely, one should, in principle, include 
deviations from this asymptotic behavior (Oxtoby and Gelbart, 1974). However, 
since the critical viscosity enhancement is weak and limited to a small range 
of temperatures anyway, we take here a practical approach by adopting (3.9) 
for all values q^>l, but treating q and (|) as adjustable parameters to be 
deduced from the experimental data. Furthermore, to connect the behavior of 
the viscosity close to the critical point with that of the normal viscosity away 
from the critical point, we neglect the critical viscosity enhancement for 
q^Sl. Thus we adopt 


T1 = n(q?)‘^ 


for 


> 1 


n = n 


for qC S 1 


(3.27) 


Using (2.20) we can also write this viscosity equation in the form 


with 


^ for X*>X* 

S = 1 for X* S X* 


<* = r/(q5 )Y/^ 


(3.28) 


(3.29) 


3.3.2 Viscosity of nitrogen 

To test the applicability of the proposed viscosity equation we 
consider a few selected fluids for which adequate experimental information is 
available. 
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Zozulya and Blagol (1974) have reported experimental data for the viscosity 
of nitrogen near the critical point. The viscosity was measured as a function 
of density directly; the usual additional uncertainties associated with the 
conversion of experimental pressures into densities are thus avoided. Zozulya 
and Blagoi determined the viscosity as a function of density along 14 isotherms 
covering the range 126.16 K£T£ 135.00 K. The critical temperature was estimated 
as T(, = (126 . 21 ± 0 . 01) K. A critical viscosity enhancement was observed in a 
temperature range of about 3 K above the critical temperature at densities between 
250 Kg/m3 and 375 kg/m^ (Pj,= 313.9 kg/m^) . Some of the experimental data are 
shown in Fig. 9. 

To represent the viscosity in the critical region we need to combine the 
equation proposed for the viscosity ratio n/O with an equation for the normal 
viscosity r|. For this purpose we adopted a representative equation proposed by 
Hanley et al. (1974a); using this equation we converted the experimental viscosity 
data into viscosity ratios ri/H- 

The viscosity ratios p/ri are plotted in Fig. 10 as a function of the 
correlation length ^ evaluated by the method described in Section 2.6. We note 
that on a double logarithmic scale the data approach a straight line for q^ > 1. 

The slope of this line yields the exponent (j) and the Intercept with the ^-axis 
the constant q“^. We thus obtain (Basu and Sengers, 1979) 



Fig. 9 Viscosity data for nitrogen in the critical region as measured by 

Zozulya and Blagol (1974) . The curves indicate the values calculated 
from the proposed viscosity equation with T^ = 126.201 K. 
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Fig. 10 Log-log plot of the viscosity ratio ri/ri as a function of the 
correlation length ^ for nitrogen in the critical region. 


q ^ = 22 X 10 c|) = 0.057 (1^= 126.201 K) , (3.30) 

so that, with the parameters in Table IV, 

X* = 10.1 , V<J)/Y = 0.0303 . (3.31) 

o 

The viscosities, calculated from (3.30) with the parameter values (3.31), are 
represented by the curves in Fig. 9. For further details the reader is 
referred to a separate publication (Basu and Sengers, 1979). 

3.3.3 Viscosity of steam 

The viscosity of steam in the critical region was investigated 
experimentally by Rivkin et al. (1975). They determined the kinematic 
viscosity r)/p as a function of temperature and pressure. We converted the 
experimental pressures into values for the density p using the equations 
presented in Section 2.5 and assuming in the analysis T^= 647.073 K as discussed 
elsewhere (Basu et al. , 1980b; Watson et al. , 1980b). The values thus deduced 
for the dynamic viscosity as a function of p are shown in Fig. 11. 
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Fig. 11 Viscosity data for steam in the critical region (Rivkin et al. 1975; 
Watson et al., 1980). The curves indicate the values calculated 
from the proposed viscosity equation with 1^=647.033 K. 

To represent the normal viscosity ri we adopted an equation analogous to 
the equation presented in Section 3.2.3 for the normal thermal conductivity 
A of steam 


with 


n = n^(T) exp 


pE E b?. d-1)^ (p- 1)^ 

i=0 j=0 


- - 1/2 ^ ^k 

n (T) = E 

° k=0 T 


(3.32) 


(3,33) 


and where T and p are the dimensionless variables earlier defined in (3.24). 
This equation has the same form as the viscosity equation adopted by the 
International Association for the Properties of Steam (Aleksandrov, 1977; 
Hendricks et al., 1977; Nagashima, 1977; Scheffler et al., 1978; Schmidt, 1979 
Meyer et al., 1980). For the coefficients a^ in (3.32) we retained the values 
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earlier determined by Aleksandrov et al. (1975) and incorporated in the current 
international viscosity equation; they are reproduced in Table VIIIA. However, 
we found that the number of nonvanishing coefficients b^^ could be reduced 
significantly as compared to the number of terms retained in the international 
viscosity equation. The coefficients bjj were redetermined in collaboration 
with Watson and they are listed in Table'^VIIIB (Watson et al., 1980a, b). 


Table VIII Coefficients in the equation for the 
normal viscosity fj of steam 


VIIIA The coefficients aP 

k 



VIIIB The coefficients b^ . 

ij 



0 

1 

2 

3 

4 

5 

0 

+0.513 2047 

+0.320 5656 

0.0 

0.0 

-0.778 2567 

+0.188 5447 

1 

+0.215 1778 

+0.731 7883 

+1.241 044 

+1.476 783 

0.0 

0.0 

2 

-0.281 8107 

-1.070 786 

-1.263 184 

0.0 

0.0 

0.0 

3 

+0.177 8064 

+0.460 5040 

+0.234 0379 

-0.492 4179 

0.0 

0.0 

4 

-0.041 7661 

0.0 

0.0 

+0.160 0435 

0.0 

0.0 

5 

0.0 

-0.015 7839 

0.0 

0. 0 

0.0 

0.0 

6 

0.0 

0.0 

0.0 

-0.003 6295 

0.0 

0.0 
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A plot of the viscosity ratios Ti/t), thus deduced from the experimental 
data, as a function of the correlation length ^ is shown in Fig. 12. 



Fig. 12 Log-log plot of the viscosity ratio ri/ri as a function of the 
correlation length ^ for steam in the critical region. 


The line corresponds to the parameters (Basu et al, 1980b) 

= 26.6 X 10“^°m, ^ = 0.05 (T^= 647.073 K) , (3.34) 

which, when combined with the parameters presented in Table V, implies 

X* = 21.4 V())/y = 0.0263 . (3.357 

The viscosities, calculated from (3.28) with the parameter values (3.35) are 
represented by the curves in Fig. 11. A comprehensive comparison of this 
viscosity equation with experimental viscosity data of water and steam over 
a large range of temperatures and densities is presented elsewhere (Watson 
et al. , 1980b) . 
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4. 


Discussion 


4.1 Recommendations for calculating the thermal conductivity in the 
critical region . 

In this report an attempt was made to formulate representative 
equations that describe the singular behavior of the transport properties in 
terms of the singular behavior of the thermodynamic properties. 

The thermal conductivity exhibits a critical enhancement in a large range 
of temperatures and densities around the critical point. In order to estimate 
this enhancement the thermal conductivity A was separated into a normal 
contribution A and a critical enhancement AA 

A = A + AA . (4.1) 


To estimate the critical enhancement we introduced an equation of the form 
(see (3.15)) 


AA 


A Vc /T*\ /gp*\ (y-V)/y 

6Trn VpV V9T*;p/ T 


P* 


exp 


[- 


{A^(AT*)^ + B^(Ap*)^ 


>]• 


(4.2) 


This equation relates the thermal conductivity enhancement at a given temperature 
and density to the pressure coefficient (8P/3T)p, the symmetrized compressibility 
Xt and the (normal) shear viscosity D- 

A calculation of the thermal conductivity enhancement AA as a function of 
temperature and pressure involves the following steps. 

a) Calculate the density at the given temperature and pressure from 
the equation of state. 

b) Calculate (3P/3T)p and Xt the given temperature and density from 
the equation of state. 

c) Calculate the (normal) viscosity n at the given temperature and density. 
In this procedure the critical viscosity enhancement may be neglected 
and most available engineering equations for the viscosity r| will be 
adequate for the purpose. 

In this report we formulated a detailed procedure for calcillating AA of 
carbon dioxide and steam in the critical region. This procedure uses available 
scaled equations of state for these substances in the critical region. 

A recommended procedure yielding estimates for the thermal conductivity 
of other fluids is presented in Table IX. We emphasize that this procedure is 
an approximate method for predicting the thermal conductivity enhancement. In 
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particular, the constants and in (4.2) will depend somewhat on the choice 
of the fluid. In the absence of further information we recommend that the 
values of these constants adopted for CO 2 be used for the other fluids as well. 
The recommended procedure was used by Hanley et al. (1974, 1975, 1977) to 
estimate the thermal conductivity enhancement for a number of fluids. 


Table IX. Procedure for estimating the critical thermal conductivity 
enhancement at given temperature and pressure 

1. Calculate p* and (9P*/9T*)p from best available equation of state. 

2. If p* and T* are outside the region bounded by (2.17), calculate 
from best available equation of state. 

If p* and T* are inside the region bounded by (2.17), calculate 
from (2.11) with parameter values from Table IV. 

3. Calculate 5^ = (a/k)'^^^ with parameter values from Table IV. 

4. Calculate the viscosity ri from the best available equation for the 
normal viscosity, neglecting any critical viscosity enhancement 
effects . 

5. Calculate AX from (4.2) with the values n=0, A=1.02, A^=39.8, 
B^=5.45 earlier adopted for CO^ (see (3.21)). 


Because of the divergent nature of the thermal conductivity at the critical 
point, the calculated thermal conductivities in the near-critical region are 
sensitive to the choice for the critical temperature and pressure. Roughly 
speaking, if the critical temperature is uncertain by a relative amount 6T*, 
then the calculated values become inaccurate in a range bounded approximately by 

I AT* I 5 |106T*|, 1Ap*| 5 |106T*|P . (4,3) 

Since at the critical point (9P*/9T*)p ~ 6, a relative error 6P* in the 
critical pressure has about the same effect as 1/6 times the same relative error 
6T* in the critical temperature. 
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4.2 Recommendations for calculating the viscosity in the critical region 


The critical viscosity enhancement is weak and can only be observed at 
temperatures within a few percent from the critical temperatures. Hence, the 
effect may be neglected in most engineering applications as was done in the 
international representation recently adopted for the viscosity of steam 
(Aleksandrov, 1977; Hendricks et al, 1977, Nagashima, 1977; Scheffler et al, 
1978; Schmidt, 1979; Meyer et al., 1980). In this report we showed how the 
critical viscosity enhancement can be accounted for in the case of nitrogen 
and steam. For other fluids it only makes sense if the density p and the 
normal viscosity r| in the critical region are known with considerable accuracy. 

An estimate of the critical viscosity enhancement can be obtained from 
(3.27) 


iln(n/n) = (|)£n(q5)0(q?-l) , (^.^) 

where 0(z) is the Heaviside step function such that 0(z) = 1 for z > 0 and 
0(z) = 0 for z-0. The correlation length ^ can be calculated from (2.20). In 
the absence of more detailed information, we recommend (f)-0.05, while q may be 
estimated from an approximate equation (Basu and Sengers, 1979) 


_ / 9P* 

‘I 8Trfi X Ut* 

c c \ 


(4.5) 


P 5, 


where n and A are the normal viscosity and normal thermal conductivity at the 
critical poinS. Using the rule of thumb (3P*/3T*) - 6, we approximate (4.5) 
by P 


q - 


9k^P r 

B c 

2irr| X 

C C O 


9k„P k 
B c 

2irri X a^^ 

c c o 


(4.6) 


The procedure is summarized in Table X. 
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Table X. Procedure for estimating the critical viscosity enhance 

ment at given temperature and density 


Calculate r| from best available equation for the normal viscosity. 

If p* and T* are outside the region bounded by (2.17), take ri=r) 
and stop. 

If p* and T* are inside the regions bounded by (2.17), continue. 

Calculate 5=?g(X^/r) C^()(^ a/k)'^^^ from (2.11) with parameter 

values from Table IV. 


Calculate and from best available equations for the normal 
thermal conductivity and viscosity. 


5. Calculate q from (4.6). 


6. Calculate ri=p(q5) ' wherever q5>l. 
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4.3 Outlook 


In this report we have presented simple approximate equations for the 
compressibility, the thermal conductivity and the viscosity of fluids in the 
critical regions. There are a number of ways available by which the quality 
of the representative equations can be Improved. 

First, it is possible to formulate a revised and extended scaled equation 
of state with a priori known theoretical values of the critical exponents and 
valid in a larger range than the scaled equation of state considered in this 
report (Balfour et al, 1977, 1978). Parameter values for such a revised and 
extended scaled equation of state have been determined in the case of steam 
(Balfour et al. 1980) and in the case of ethylene (Hastings et al., 1980) and 
parameter values in the case of other fluids are expected to become available 
in the future. 

For a completely satisfactory description one would like to formulate a 
global equation of state that smoothly connects the scaling laws in the 
critical region with the known analytic surface valid away from the critical 
region. This problem requires further study as discussed elsewhere (Levelt 
Sengers and Sengers, 1980). 

In order to connect the critical enhancement in the transport properties 
with the normal behavior of the transport properties we introduced an 
empirical crossover function (3.13) for the thermal conductivity and a 
Heaviside step function (see (4.4)) for the viscosity. It is possible to 
Introduce better crossover functions from a more detailed examination of 
the mode-coupling equation for the critical enhancement of the thermal 
conductivity and the viscosity. This work is currently in progress."^ 


t 


This work is done in collaboration with Bhattachar j ee and Ferrell. 


42 



References 


Aharony, A. and G. Ahlers (1980), Phys. Rev. Lett. 782. 

Aleksandrov, A. A. (1977), Thermal Engineering 2 ^ (4), 76. 

Aleksandrov, A. A., A. I. Ivanov and A. B. Matveev (1975), Thermal 
Engineering 22^ (4), 77. 

Aleksandrov, A. A. and A. B. Matveev (1976), "Equation for Thermal Conductivity 
of Water in the Wide Range of Temperatures and Pressures", Report submitted to 
the Special Committee of TAPS (International Association for the Properties of 
Steam) . 

Baker, G. A., B. G. Nickel and D. I. Meiron (1978), Phys. Rev. B 2 ^, 1365. 

Balfour, F. W. , J. V. Sengers, M. R. Moldover and J. M. H. Levelt Sengers (1977), 
in Proceedings 7th Symposium on Thermophysical Properties , A. Cezairliyan, ed . 
(Americal Society of Mechanical Engineers, New York), p. 786. 

Balfour, F. W. , J. V. Sengers, M. R. Moldover and J. M. H. Levelt Sengers (1978), 
Phys. Lett. 65A , 223. 


Balfour, F. W. , J. V. Sengers and J. M. H. Levelt Sengers (1980), in Water and 
Steam: Their Properties and Current Industrial Applications , J. Straub and K. 
Scheffler, eds. (Pergamon Press, Oxford and New York), p. 128. 

Barber, M. N. (1977), Phys. Reports 29C , 1. 


Basu, R. S. and J. V. Sengers (1977), in Proc . 7th Symposium on Thermophysical 
Properties, A. Cezairliyan, ed. (American Society of Mechanical Engineers, 

New York), p. 822. 

Basu, R. S. and J. V. Sengers (1979), J. Heat Transfer, Trans ASME 101 , 3; 

101 , 575. 


Basu, R. S. and J. V. Sengers (1980) in Thermal Conductivity - 16 , D. C. 

Larsen, ed. (Plenum Publ. Corp., New York), in press. 

Basu, R. S., J. V. Sengers and J. T. R. Watson (1980), Intern. J. Thermophys. 

1, 33, 

Becker, H. and U. Grlgull (1975), Maschinenmarkt 50. 

Becker, H. and U. Grigull (1978) , Warme - und Stof fubertragung lA, 9. 

Benedek, G. B. (1968), in Statistical Physics, Phase Transitions and Superfluidity , 
Vol. 2, M. Chretien, E. P. Gross and S. Deser, eds. (Gordon and Breach, New 
York) , p. 1, 


43 


References (cont’d.) 

Burstyn, H. C., J. V. Sengers and P. Esfandlari (1980), Phys. Rev. k 282. 
Calmettes, P. (1977), Phys. Rev. Lett. ^9^, 1151. 

Cummins, H. Z. and H. L. Swinney (1970) in Progress In Optics , Vol. 8, E. Wolf, 
ed. (North Holland Publ. Comp., Amsterdam), p. 135. 

De Domlnicls, C. and L. Peliti (1978), Phys. Rev. B IB, 353. 

De Pasquale, F. , P. Tartaglia and P. Tombesi, (1978), J. Phys. A 13^, 2033. 

Dlller, D. E. , Hanley, H. J. M. and Roder, H. M. (1970), Cryogenics 1^, 268. 

Domb, C. and M. S. Green, eds. (1976), P hase Transitions and Critical Phenomena , 
Vol. 6 (Academic Press, New York). 

Douglass, T. B. (1969), J. Res. Natl. Bur. Stand. 73A , 541. 

Ernst, M. H. (1978), in Fundamental Problems in Statistical Mechanics IV , E. G. D. 
Cohen and W. FisEdon, eds. (Zaklad Narodowy um Ossolinskich-Wydawnictwo, 

Wroclaw, Poland), p. 161. 

Ferrell, R. A. (1970), Phys. Rev. Lett. 24 , 1169. 

Fleury, P. A. and J. P. Boon, (1973) Adv. Chem. Phys. 2^, 1. 

Greer, S. C. (1978), Acc. Chem. Res. 11, 427. 

Griffiths, R. B. (1967), Phys. Rev. 158 . 176. 

Guttinger, H. and D. S. Cannell (1980), Phys. Rev. A. 285. 

Gunton, J. D. (1979), in Dynamical Critical Phenomena and Related Topics , C. P. 
Em, ed. (Springer Verlag, Berlin, Heidelberg, New York), p. 1. 

Gunton, J. D. and K. Kawasaki, (1976), Progr. Theor. Phys. 61. 

Hanley, H. J. M. , R. D. McCarty and E. G. D. Cohen (1972), Physica 6^, 322. 

Hanley, H. J. M. , R. D. McCarty and W. M. Haynes (1974a), J. Phys. Chem. Ref. 

Data 979. 

Hanley, J. H. M. , K. E. Gubbins and S. Murad (1977), J. Phys. Chem. Ref. Data 
6, 1167. 

Hanley, H. J. M. , R. D. McCarty and W. M. Haynes (1975), Cryogenics 1^, 413. 


44 



References (cont’d.) 


Hanley, H. J. M. , R. D. McCarty and J. V. Sengers (1974b), "Viscosity and 
Thermal Conductivity Coefficients of Gaseous and Liquid Oxygen", NASA 
Contractor Report 2440 (National Aeronautics and Space Administration, 
Washington , D . C . ) • 

Hanley, H. J. M. , J. V. Sengers and J. F. Ely (1976), in Thermal Conductivity - 
14 . P. G. Klemens and T. K. Chu, Eds. (Plenum Publ. Corp., New York), p. 383. 

Hastings, J. R., J. M. H, Levelt Sengers and F. W. Balfour (1980), J. Chem. 
Thermodynamics 3^, 1009. 


Hendricks, R. C., R. B. McClintock and G. J. Silvestri (1977), J. Eng. Power, 
Trans. ASME 664. 

Hohenberg, P. C. and M. Barmatz (1972), Phys. Rev. A 289. 

Hohenberg, P. C. and B. I. Helper in (1977), Rev. Mod. Phys. 4_9, 435. 

Hubbard, J. and P. Schofield (1972), Phys. Lett. 40A , 245. 

IFC (1968) , "The 1968 IFC Formulation for Scientific and General Use" (American 
Society of Mechanical Engineers, New York). 

Kadanoff, L. P. and J. Swift (1968), Phys. Rev. 166 , 89. 

Kawasaki, K. (1970), Annals of Physics^, 1. 

Kawasaki, K. (1971), in Critical Phenomena , Proc. Intern. School "Enrico 
Fermi", Course LI, M. S. Green, ed. (Academic Press, New York), p. 342. 

Kawasaki, K. (1976) in Phase Transitions and Critical Phenomena , Vol. 5A, 

C. Bomb and M. S. Green, eds. (Academic Press, New York), p. 165. 

Kawasaki, K. and J. D. Gunton (1976), Phys. Rev. B 4658. 

Kestin, J. (1978), Mechanical Engineering 100 (8), 46. 

Kestin, J. (1980) , in Wa ter and S team; Their Properties and Current Industrial 
Applications . J. Straub and K. Scheffler, eds. (Pergamon Press, Oxford and New 
York), p. 283. 

Le Guillou, J- C. and J. Zinn-Justin (1977), Phys. Rev. Lett. 95. 

Le Guillou, J. C. and J. Zinn-Justln (1980), Phys. Rev. B 21 , 3976. 

Le Neindre, B., R. Tufeu, P. Bury and J. V. Sengers (1973), Ber. Bunsenges. 
physik. Chemie Tl_, 262. 

Levelt Sengers, J. M. H. (1974), Phys lea 73^, 73. 


45 


L 


Levelt Sengers, J. M. H. (1977a), in Proc. 7th Symposium on Thermophysical 
Properties, A. Cezairliyan, ed. (American Society of Mechanical Engineers, 
New York) , p. 766. 

Levelt Sengers, J. M. H. (1977b), in Proc. 7th Symposium on Thermophysical 
Properties, A. Cezairliyan, ed. (American Society of Mechanical Engineers, 
New York) , p. 774. 

Levelt Sengers, J. M. H. and W. T. Chen (1972), J. Chem. Phys. 56 , 595. 

Levelt Sengers, J. M. H., W. L. Greer and J. V. Sengers (1976), J. Phys. 
Chem. Ref. Data _5, 1. 


Levelt Sengers, A., R. Hocken and J. V. Sengers (1977), Physics Today 3^0 (12), 42. 

Levelt Sengers, J. M. H. and J. V. Sengers (1975), Phys. Rev. A 1^, 2622. 

Levelt Sengers, J. M. H. and J. V. Sengers (1980), in Perspectives in 
Statistical Physics , M. S. Green Memorial Issue of "Studies in Statistical 
Mechanics", H. J. Raveche, ed. (North-Holland Publ. Comp., Amsterdam), in 
press . 

Ley-Koo, M. and M. S. Green (1977), Phys. Rev. A 1^, 2483. 

Mazenko, G. F. (1978), in Correlation Functions and Quasiparticle Interactions , 

J. W. Halley, ed. (Plenum Press, New York), p. 1. 

McIntyre, D. and J. V. Sengers (1968), in Physics of Simple Liquids , H. N. V. 
Temperley, J. S. Rowlinson and G. S. Rushbrooke, eds. (North Holland Publ. 

Comp., Amsterdam), p. 447. 


Meyer, C. A., R. B. McClintock, G. J. Silvestrl and R. C. Spencer (1980), ASME 
Steam Tables , 4th ed. (American Society of Mechanical Engineers, New York). 

Michels, A., B. Blaisse and C. Michels (1937), Proc. Roy. Soc. (London) A 160 , 
358. 


Michels, 

A. 

and 

Michels , 

A. 

and 

Michels , 

A. 

and 

Michels , 

A. 

9 J • 

Michels , 

A. 

and 

Moldover 

, M 

. R. 


S. R. de Groot (1948), Appl. Sci. Res. A 94. 

C. Michels (1935), Proc. Roy. Soc. (London ) A 153 , 201. 

C. Michels (1937), Proc. Roy. Soc. (London) A 160 , 348. 

V. Sengers and P. S. van der Gulik (1962), Physica 28 , 1216. 
J. C. Strijland (1952), Physica 18, 613. 

(1974), J. Chem. Phys. ^1^, 1766. 


46 



Holdover, M. R. (1978), J. Res. Natl. Bur. Stand. 329. 

Mountain, R. D. (1966), Rev. Mod. Phys. 38, 205. 

Murphy, T. A., J. V. Sengers and J. M. H. Levelt Sengers (1973), in Proc . 

6th Symposium on Thermophysical Properties , P. E. Liley, ed. (American 
Society of Mechanical Engineers, New York), p. 180. 

Murphy, T. A., J. V. Sengers and J. M. H. Levelt Sengers (1975), in Proc . 

8th Int ernational Conference on the Properties of Water and Steam , P. Bury, 
H. Perdon and B. Vodar, eds. (Editions Europeennes Thermiques et Industries, 
Paris) , p . 603 . 

Nagashima, A. (1977), J. Phys. Chem. Ref. Data 1133; J_, 1755. 

Nlcoll, J. F. (1980), Phys. Lett. 76A , 112. 


Nicoll, J. F. and T. S. Chang (1978), Phys. Rev. A 2083. 
Ohta, T. (1975), Progr. Theor. Phys. 1566. 

Ohta, T. (1977), J. Phys. C 791. 

Ohta, T. and K. Kawasaki (1976), Progr. Theor. Phys. 5^, 1348. 
Oxtoby, D. W. and W. M. Gelbart (1974), J. Chem. Phys. 61^, 2957. 
Perl, R. and R, A. Ferrell (1972), Phys. Rev. A 2358. 


Pfeuty, P. and G. Toulouse (1977), Introduction to the Renormalization 
Group and to Critical Phenomena (Wiley, New York). 

Pomeau, Y. and P. Resibois (1975), Phys. Reports 19C , 63. 


Rivkin, 

S. 

L. 

and 

T. 

S . Akhundov 

(1962), 

Teploenerget ika 

9_ 

(1), 57. 

Rivkin , 

S. 

L. 

and 

T. 

S . Akhundov 

(1963), 

Teploenergetika 

10 

(9), 66. 

Rivkin, 

S. 

L. 

. T. 

S. 

Akhundov , E . 

A. Kremenevskaya and N, 

, N 

. Assadulaieva 


(1966), Thermal Engineering 3^ (4), 77. 

Rivkin, S. L. , A. A. Aleksandrov and E. A. Kremenevskaya (1978), Thermodynami c 
Derivatives for Water and Steam (Wiley, New York) . 

Rivkin, S. L. , A. Ya. Levin, L. B. Izrailevsky and K. G. Kharitonov (1975), in 
Proc. 8th International Conference on the Properties of Water and Steam , P. Bury, 
H. Perdon and B. Vodar eds. (Editions Europeennes Thermiques et Industries, 
Paris), p, 153. 


47 




Scheffler, K. , N. Rosner, J. Straub and U. Grigull (1978), Brennst - Warme - 
Kraft 3G, 73. 

Scheffler, K. , N. Rosner, J. Straub and U. Grigull (1979), Brennst - Warme - 
Kraft 326. 

Schmidt, E. (1979), Properties of Water and Steam in SI - Units , 2nd ed., 

U. Grigull, ed. (Springer Verlag, Berlin, Heidelberg, New York). 

Schmitt, H. and H. Beer (1977), Brennst - Warme - Kraft 15. 

Schofield, P. (1969), Phys. Rev. Lett. 2^, 606. 

Schofield, P., J. D. Litster and J. T. Ho (1969), Phys. Rev. Lett 2^, 1098. 

Scott, R. L. (1978), in Chemical Thermodynamics , Vol. 2, Specialist Periodical 
Report, M. L. McGlashan, ed . (Chemical Society, London), p. 238. 

Sengers, J. V. (1966), in Critical Phenomena , NBS Misc. Publ . 273, M. S. 

Green and J. V, Sengers, eds. (U.S. Gov't. Printing Office, Washington, D.C.), 
p. 165. 

Sengers, J. V. (1968), in Recent Advances in Engineering Science , Vol. 3, 

A. C. Eringen, ed . (Gordon and Breach, New York, 1968), p. 153. 

Sengers, J. V. (1971), in Critical Phenomena, Proc. Intern. School "Enrico 
Fermi" , Course LI, M. S. Green, ed . (Academic Press, New York), p. 445. 

Sengers, J. V. (1972), "Transport Properties of Gases and Binary Liquids near 
the Critical Point", NASA Contractor Report 2112 (National Aeronautics and 
Space Administration, Washington, D.C.). See also: Ber . Bunsenges. physik. 

Chemie 7^, 234 (1972) . 

Sengers, J. V. (1973), in Transport Phenomena - 1973 , AIP Conference Proceedings 
No. 11, J. Kestin, ed. (American Institute of Physics, New York), p. 229. 

Sengers, J. V. and R. S. Basu (1977), "The thermal conductivity of steam in 
the zero density limit as a function of temperature". Technical Report BN 852 
(Institute for Physical Science and Technology, University of Maryland, 

College Park, MD) . 

Sengers, J. V., R. S. Basu, J. T. R. Watson and R. C. Hendricks (1981), J. Phys. 
Chem. Ref. Data, to be published. 

Sengers, J. V. and P. H. Keyes (1971), Phys. Rev. Lett. 2 §_, 70. 

Sengers, J. V. and J. M. H. Levelt Sengers, (1977), "Concepts and methods for 
describing critical phenomena in fluids", NASA Contractor Report 149665 (NASA 
Lewis Research Center, Cleveland, OH). With some minor modifications this 
report was published in Progress in Liquid Physics , C. A. Croxton, ed. (Wiley, 
New York, 1978), p. 103. 


48 



Sengers, J. V. and M. R. Holdover (1978), Phys. Lett. 66A , 44. 

Slggia, E. D. , B. I. Halperin and P. C. Hohenberg (1976), Phys. Rev. B 13 , 

2110. 

Sirota, A. M. , V. I. Latunln and G. M. Belyaeva (1973), Thermal Engineering 

20 (8), 7. 

Sirota, A. M. , V. I. Latunin and G. M. Belyaeva (1974), Thermal Engineering 

21 (10), 70. 

Sirota, A. M. , V. I. Latunin and G. M. Belyaeva (1976), Thermal Engineering 
^ (1), 72. 

Sirota, A. M. , V. I. Latunin, G. H. Belyaeva and I. I. Bol'dshtein, (1976), 
Thermal Engineering 23^ (6), 75. 

Vlcentini-Missoni, M. , J. M. H. Levelt Sengers and M. S. Green, (1969a), J. 

Res. Natl. Bur. Stand. 73A , 563. 

Vicentlni-Missonl, M. , J. M. H. Levelt Sengers and M. S. Green (1969b), Phys. 

Rev. Lett. 2^, 389. 

Wallace, D. J. and R. K. P. Zia (1978), Rep. Prog. Phys. 4T, 1. 

Watson, J. T. R. (1976a), "Dynamic viscosity of carbon dioxide gas and liquid", 
ESDU Data Item No. 76021 (Engineering Sciences Data Unit, London). 

Watson, J. T. R. (1976b), "Thermal conductivity of carbon dioxide gas and 
liquid", ESDU Data Item No. 76030 (Engineering Sciences Data Unit, London). 

Watson, J. T. R. (1978), "Thermal conductivity of ice, water and steam", 

ESDU Data Item No. 78039 (Engineering Sciences Data Unit, London). 

Watson, J. T. R. , R. S. Basu and J. V. Sengers (1980a), in Water and Steam : 

Their Properties and Current Industrial Applications , J. Straub and K. Scheffler, 
eds. (Pergamon Press, Oxford and New York), p. 336. 

Watson, J. T. R. , R. S. Basu and J. V. Sengers (1980b), J. Phys. Chem. Ref. 

Data, 9_, 239. 

Waxman, M. , H. A. Davis, J. M. H. Levelt Sengers and M. Klein (1978), "The 
equation of state of Isobutane: An interim assessment", Nat'l. Bur. Stand. 

Internal Report, prepared for the Geothermal Energy Division of the U.S. 
Department of Energy. 

Wilson, K. G. and J. Kogut (1974), Phys. Reports 12C , 75. 

Zozulya, V. N. and Yu P. Blagoi (1974), JETP 39, 99. 


49 


Appendix 


Computer program for calculating the linear model parameters r and 0 


C iUEPCUTINE CPPfiK 

C ftiJTKGR T.fi.nuPPHV 

C 

c PIJPPC5E 

C EVRI.IJPTE LI HEAP KCCEL PRPRKETERS R RUE TKETR 

C FPDM GIVEN TENPEPfiTUPES RUB DENSITY 

C 

C LISRGE 

C CRLL EPPRN<ERS.ETS.NCBE.XCS.E^S-EETRS*R.TK) 

C 

C EESCPIPTICN CP PRPRKETERS 

C EPS = PKC^PHCC-t=PEEL:CEB DENSITY- 1 

C DTS = T^TC-1=PEDIJCED TENPERRTUPES- 1 

C PSfTKS - LINERP KCEEL PRRRNETERS R. Th'ETR 

C EETRS - CPITICRL EXPONENT 

C XCS.E2S - FIXED PRPRNETERS IN LINERR NOEEL 

C MODE - RN INEICRTCP 

C = 1 IF ERTR POINT IS INSIDE CXC 

C = a IF NEWTONS NETKOB CRNNDT FINE THETR 

r = 3 IF PRRRNETEPS P RNE TKETR RPE SUCCESSFULLY EVRLURTEE 

C = 4 IF TKETR WRS FOUND BUT R FRILS CONSISTENCY TEST 

C = S IF TKETR WRS FOUND BUT R BECCNES NECRTIVE 

C 

C CGNNENTS 

C TKE PROSRRK USES BCUBLE PRECISION RRITKNETIC 

C ERS.ETS. EETRS. X0S.B2S SKCULE BE PRCVIEEB BY TKE CRLLING ROUTINE 

C TKE VRLUE OF R RNE TKETR IS RETURNEE TC TKE CRLLING ROUTINE 

C TKE INPUT PRRRNRTERS SKCULE BE IN SINGLE PRECISICN 

C TKE PRCFRRK CONVERTS TKEN TO BCUBLE PRECISION RNE USES 

C DOUBLE PRECISION RRITKNETIC . RETURNS R. TKETR IN DOUBLE 

C RRECISION 

C 

C SUBROUTINES RNE SUBPR06RRNS REGUIREE 

C NONE 

C 

SUBRDUT I NE DRRRN .'DRS . DTS . NODE. XOS . BSS . BETRS . R . TK> 

INPLICIT RERL«8<R-H.C-2> 

PERL DRS.DTS.XCS.BETRS.B2S 
I>1 
C 

C INRX = ICC GIVES TKE NRXINUN NUNBER OF ITERATIONS DONE 

C ERNRX SETS TKE CRITERIR FDR CONVERGENCE 

C 

INflX-lCC 
DR « DRS 
DT = DTS 
XC - XOS 
BETR - BETRS 
B2 » E2S 
ERNRX- l.E- 13 
B-SGRT<B2> 

SGN=1. 

IFCEP.LT.C.T SGN--1. 

IFCDR.EC. C.T GO TO 2CC 
IF<ET.EC. C.> GO TO ICC 
X=DT''flESa:R>«« fl. --BETRT 
RX-RBS <X') 

IF<X.LT. (-XD1> 'GD TC 4CC 
IF ^DT.GT.C.l GO TC 5C 
C 

C NEWTONS NETKCD OF SOLUTION FOP VRRIDUS VRLUES OF ER RND DT 

C 

THETR=a.*l./B>^2. 

25 FX=RX* <E2- 1 .>*THETR*« TI.xBETRT ^XC* 1 . -E2»THETR«*2 

FXP=flX« <E2- 1 . ) «THETR4* a . --BETR- I , > x <XD«BETfl> -2. *B2»THETR 
TKETR=TKETR-FX''FXP 

IF <RBS <FX'' <FXP«THETR" .LT. ERNRX') 60 TO 3CC 
IFd.GT. INRX) 60 TO 5CC 
i-r-M 
GO TO 25 
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5C THETR 

5^ FX=ftX« <B3- 1 . > «TKETR*« <1 . /'BETR> ✓XC- 1 . +BS«THETR«*2 

FXP=RX« (Ba- 1 . ) ♦TRETR** <1 . ^BETR- 1 . > ^ OiC^BETR') +a. «Ba« THETR 
TKETR=THETR- FX^FXP 

IF(RES<FXy-(FXP«THETR>>.LT.ERKRX> GC TC 3CC 

I»I41 

IFa.GT.IKflX) 6C TC 5CC 
6C TC 55 

ICC TRETR=1./'B*S6N 
TH=TKETR 
MCEE*3 

R= f BR^THETR > ♦ ♦ < 1 . /'BETfl'J ♦XC/' <B8- 1 - > 

RETURN 
aCC TRETR-C. 

TH=THETfl 

NCEE«3 

R=DT 

IF<CT.LT.C.) 6C TC 4CC 
RETURN 

3CC TRETR=THETfl«S6N 
THi>THETfl 

DDUBR=DT ^ < 1 . -B£«THETR*«a> 

R=I!CUBR 

IFfCCUER.LT. C.) 6C TC 6CC 
DRC :RLC= <ECUER« f Ba- 1 . > /XC> ♦♦BETR*THETR 
IF fRESa.-CRCRLC^ER>,6T. I.E-IO GC TC 35C 
NCEE=3 
RETURN 
35C KCEE»4 

WRITE <6.9CCl>CRCflLC.CRtBT,R»TH»XC.BETfl.Ba, I.KCEE 
ERNftX=ERNfiX^l CCCC 

9CC1 FCRNRT^-^CPRRNTER IN TRCUBLE- NRY REQUIRE KILLING RUN',-^. 
14Ea4. 18. -'4Ca4. 18.^.ai3> 

THETR=RBS<THETfl) 

1 = 1 

IFCET.GT. C.l GC TC 55 
GC TC as 
4CC NCEE-1 
TH=S6N 

BCUBR=ET^<1.-Ba> 

R=ECUBR 
RETURN 
CC NCEE=a 

SC URITE<G.'?CCa>CR.ET.R.TH.XC.BETR.Ba. I.NCEE 
9CCa FCRKRT (' CPRRKTR FfilLURE' . 7E1 4. 8. z'. a I S') 

RETURN 
6CC NCEE*5 

GC TC 55 C 
ENE 
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